Introduction: Safety and sustainability of infrastructures which were placed in or on rock mass mainly control by geometrically size distribution and physical and mechanical characteristics of rock blocks that is created by intersection of discontinuities. hence identification of rock blocks has a key role in mechanical analysis and hydraulic behaviour of jointed rock mass. Detection process of blocks have many applications in rock mechanic which could be referred to their use in the numerical methods like discrete element method or in analysis of continuous deformation of discontinuities. As pioneer researchers, Goodman and Shi, Warburton and Heliot could be known as leaders in the field of diagnosis of rock mass blocks. Warburton provides a method based on geometric parameters of rock mass and developed a software based on it. Warburton in his work assumed discontinuities as parallel and infinite. In the earlier works, discontinuities were considered as infinite panes. So, just convex blocks were distinguishable. Concave blocks were diagnosis in more detailed researches that is created by finite discontinuities. Basically, methods based on finite planes was classified into two branches. Aforementioned branches were based on blocks detection based on topology concepts and assemble of block elements. Lin at al. presented detection method that assumed discontinuities as finite planes and worked based on topology theory. This method could realize convex and concave blocks of rock mass. Ikegawa and Hudson, Jing presented the similar methods using more accurate process. Sharma et al. presented an equation for calculating the volume of rock blocks in their work. Ferreira provided a method based on graph theory which is better than other method considering time and complicity. Based on this method, firstly vertices were detected in two dimensions and then created a graph based in vertices and edges which in next step constitute polygons that are form in two-dimension blocks. In the present research, it is developed high-speed algorithms to identify the blocks. New method was developed in MATLAB software that by assuming infinite discontinuities and inclusion of a set of joints. we have identified created blocks and calculated their volume and at last block volume histogram were draw that paves the way to obtain their distribution function. Material and methods: Infinite planes are used to simulate of discontinuities. in this study, each discontinuity is represented by a plane in a three-dimensional Euclidean space. To identify the block, a certain volume of rock mass space should be considered as study region. The studied volume is called domain. By the intersection of discontinuity planes in space, rocky blocks are created in the domain. First, vertices should be recognized at first as first step in block detection. Then, edges are diagnoses and after that it's time to specify the polygons and finally, polyhedron or blocks are obtained by joining edges together. Each vertex in space is created by the intersection of three nonparallel planes. In fact, the vertex is the interface of three planes in the Euclidean space. The next element in the block metric process is the diagnosis of the edges or the blocks' edges. All edges are sections on the lines which created by the intersection of the planes in space. first the parallel vector of all the lines resulting from the intersection of the pair of planes is obtained. After detection of edges, it’ s time to identify polygons that form key element of blocks. Each polygon of a block is formed from their constituent unit. In this step, polygons belong to each discontinuity plane is identified separately. Some edges are determined that are start from the end of selected edge between other edges. In this state, if there is just one edge, that edge is record as the next edge of first polygon. If there is more than one edge from the edge of the selected edges, the angle is calculated between each possible of end edge with the selected edge. In the next step, it’ s time to diagnosis polyhedrons that have created by discontinuities intersection. In the previous step, possible polygons were obtained for each discontinuity. In this stage, it is used the principle which is designed this algorithm that two polygons that formed a block have a common edge. So, the first polygon of first discontinuity is consider as first polygon of first block to recognize block. Results and discussion: According to the developed algorithm, MATLAB software was used to model the discontinuities. The computational and graphic capabilities of this software have created a lot of attractions for most researchers to use its potential. The strengths of this software are high computing power with its graphical accuracy. The code developed in MATLAB is called RockBlock2 that is designed using a graphical user interface (GUI) to make it easy to use. To illustrate how the program works, there are 29 discontinuities given to the program. The program first takes the dip and dip direction of discontinuities along with the desired point on it and calculates the parameters that make up the equation of discontinuity planes. Input data is stored in a separate Excel file that was previously introduced to the program. In the next step, the program attempts to identify the vertices. The program stores the coordinates of each corner, with the assignment of a number to it, in the matrix of the corners, which is in fact the Excel file that was previously introduced to the program to use in the next steps, after recognizing vertices on the area. Identifying the edges is the next step that the program done. At this stage, the program begins to identify each single edge using the data from the previous step that means the coordinates of the corners and the algorithm defined. The coordinates of the beginning and end of each edge along with its number are stored and maintained in the edge matrix in the Excel file format. In the stage of identifying the polygons, the polygons are formed by joining the edges together. This matrix is a special matrix that its matrix matrices are matrix itself. The matrix of polygons is a row matrix; whose number is the number of discontinuities. Because, as it mentioned in the chapter of the algorithm, the polygons are found by separation of discontinuities. Therefore, each column of the polygons matrix is consisting of faces that are on a certain discontinuity. The next step begins the process of identifying the blocks, or the same polygons by the program. At this step, the program starts the identification process using the features found in the previous step and the algorithm defined for it. At this stage, the identified blocks are stored in the blocks matrix. By identifying blocks, the program calculates the volume of each block and finally draw its volume histogram. In fact, a volume histogram is presented to illustrate how the block volume is distributed. Obtaining the distribution of blocks or, in other words, achieving a block probability distribution function is an essential step in the behavior of rock mass. Because one of the most important consequences of the presence of discontinuities is the fragmentation of the rock material under the block intervals. By having the block distribution function, it is possible to produce a blockbuster method using random methods, such as Monte Carlo, and to analyze it in various and arbitrary modes. Conclusion: To identify and study the rocky blocks created by discontinuities, a hierarchical algorithm was designed and developed in MATLAB software. This algorithm identifies and records blocks, consisting of blocks, edges, and facets of the blocks forming components, including stone blocks. This algorithm, which is written for user-friendly ease with the use of graphical coding capabilities, shows a very fast performance using the parallel computing power of MATLAB software. The developed code calculates the dip and dip direction of discontinuities using the geometric properties, and calculates the blocks created in three dimensions and calculates their volume. This histogram code displays the calculated volumes. The results show that the developed code with its fast performance, while identifying the blocks, calculates and records their volumes without errors. The ability to display the step-by-step process of identifying blocks is one of the clear features of this code. Information about edge is also records and is available for auxiliary applications. Histogram of block volume is one of the most important results of the developed code, which can have different applications. Identification of created rocky blocks is used both in the stability analysis and rock mass simulations such as Discrete Fracture Network modeling. Determination of block volume distribution function which is done using histogram is one of the most important uncertainties in three-dimensional rock masses behavior that can play a key role in optimizing the design of structures involved in rock mass. Therefore, considering the key role of blocks volume, identifying and calculating block volumes and, consequently, plotting their histogram and determining the distribution function governing them, has a key role in the static and dynamic analysis of rock base structures.